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SUMMARY 


This  report  pertains  to  the  overall  problem  of  estimating  damping  due  to 
dry  friction  forces  induced  at  interfaces  of  vibrating  components.  The  phase 
of  the  effort  reported  here  contains  the  summary  of  a  literature  survey  and 
identifies  the  scope  of  research  needed. 


The  survey  clearly  indicated  the  complex  processes  involved  during  rub¬ 
bing  of  one  component  relative  to  another  bringing  into  play  a  host  of  param¬ 
eters.  Very  limited  data,  especially  of  a  fundamental  nature,  appears  to  be 
available  in  regard  to  dry  friction  in  a  vibratory  environment. 


An  analytical  and  corresponding  experimental  plan  is  presented  and  dis¬ 
cussed.  These  plans  include  (a)  considerations  of  new  laws  of  friction  appro¬ 
priate  to  vibration  analyses  and  (b)  corresponding  basic  experiments  from 


which  a  data  base  can  be  developed  for  use  in  design  of  structural  compo¬ 
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The  principal  contents  of  this  report  were  'presented  in  the  form  of  a 
paper  entitled  "Characterist ics  of  Dry  Friction  Damping",  at  the  Air  Force 
sponsored  Vibration  Damping  Workshop,  held  in  Long  Beach,  California, 
February  27-29,  1984. 


This  annual  report  covers  the  period  between  March  1,  1983  and 
March  1,  1984.  The  work  reported  herein  was  under  sponsorship  of  the  AFOSR 
under  the  technical  direction  of  Dr.  A.  Amos. 
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1 .  INTRODUCTION 


The  structural  integrity  of  engineering  systems  depends  to  a  large  extent 
on  the  levels  of  energy  dissipation  that  can  occur  when  the  system  is 
subjected  to  time  dependent  forces.  One  of  the  most  important  sources  of 
energy  dissipation  in  built-up  structures  pertains  to  dry  friction  at  joints 
and  at  interfaces  in  contact  with  each  other.  When  materials  with  very  low 
inherent  damping  are  used  in  environments  offering  essentially  no  aerodynamic 
damping  (such  as,  for  example,  in  space  structural  systems),  the  ability  of 
the  structures  to  withstand  resonant  vibration  depends  almost  entirely  on  the 
extent  of  friction  damping.  Even  in  the  design  of  components  which  can 
generate  certain  levels  of  aerodynamic  damping,  it  is  essential  to  be  able  to 
estimate  accurately  the  levels  of  non-aerodynamic  sources  of  damping  because 
unlike  the  latter,  the  former  can  be  either  positive  or  negative.  For 
example,  the  aeroelastic  instabilities  of  aircraft  structures  such  as  airplane 
wings  or  turbomachinery  blades  are  attributed  to  negative  aerodynamic  damping 
developing  in  the  system  as  a  result  of  its  interaction  with  the  air  forces 
causing  the  system  to  vibrate.  In  such  situations,  escalation  of  vibratory 
stresses  to  dangerous  levels  is  prevented  only  through  the  contribution  to 
damping  from  sources  such  as  friction.  In  some  instances,  artificial  devices 
are  designed  and  introduced  into  the  structural  system  with  the  sole  purpose 
of  developing  friction  forces  leading  to  damping  when  the  dynamics  of  the 
component  demands  it.  Thus,  there  is  a  clear  need  to  enhance  the  level  of 
understanding  in  the  broad  area  of  dynamics  of  surfaces  in  contact. 

Friction  between  contacting  interfaces  which  undergo  relative  vibratory 
motion  is  known  to  dissipate  the  energy  of  vibration  resulting  in  damped 
oscillations.  The  phenomenon  of  friction  between  contacting  surfaces  is 
probably  the  most  elusive  physical  mechanism  that  defies  clear  comprehension. 
The  complexity  of  the  phenomenon  of  friction  damping  arises  from  the  varia¬ 
tions  in  the  type  of  time  dependent  motions  developing  at  an  interface.  These 
variations  in  the  relative  motion  at  the  contacting  surfaces  span  the  extremes 
between  microslip  and  gross  motion  and  include  local  slip,  stick-slip  motion, 
chatter,  etc.  The  parameters  that  control  the  resulting  motion  include  the 
normal  forces  holding  the  surfaces  together,  their  distribution,  properties  of 
materials  in  contact,  surface  treatments,  temperature,  frequency  of  vibration, 
level  of  vacuum,  coefficient  of  dynamic  friction,  etc.  The  phenomenon  is 
clearly  nonlinear  and  the  feasibility  of  linearization  needs  to  be  estab¬ 
lished  . 

In  the  context  of  vibration  engineering,  the  basic  requirements  are  to 
(a)  quantify  the  nature  and  magnitude  of  friction  forces  that  are  manifest 
between  contacting  interfaces  of  vibrating  co,.  oonents  ,  (b)  quantify  the  nature 
and  magnitude  of  vibratory  motion  at  these  interfaces,  and  (c)  predict  the 


extent  of  damping  that  may  be  present.  Within  this  context,  friction  forces 
are  considered  to  be  useful,  i.e.  they  control  vibratory  amplitudes  which 
otherwise  may  escalate.  On  the  other  hand,  any  consideration  of  friction 
forces  cannot  ignore  the  influence  of  these  forces  on  wear  of  the  components 
resulting  in  loss  of  useful  life  of  machines.  The  phenomena  of  friction  and 
wear  are  thus  inseparable.  The  emphasis  in  this  study,  however,  is  in  out¬ 
lining  those  aspects  of  friction  forces  that  pertain  to  vibration  damping. 


2.  REVIEW  OF  THE  LITERATURE 


In  order  to  assess  the  state-of-the-art  in  the  area  of  damping  due  to  dry 
friction,  a  literature  search  was  initiated.  It  became  clear  that  the  subject 
matter  could  be  classified  into  at  least  three  specialized  disciplines 
comprising  efforts  in  the  areas  of  (a)  surface  science  (b)  characterization  of 
friction  forces  at  mating  surfaces  and  (c)  damping  technology.  Further,  there 
appears  to  be  very  limited  data,  especially  of  a  fundamental  nature,  obtained 
in  a  vibratory  environment  in  which  friction  forces  play  a  significant  role. 
The  developments  in  each  of  these  areas  have  occurred  somewhat  independent  of 
each  other  because  each  represents  an  area  important  in  itself  providing 
considerable  research  opportunities.  Future  developments  must,  however, 
attempt  to  integrate  these  efforts  in  order  to  be  able  to  obtain  accurate  and 
reliable  estimates  of  energy  dissipation  in  vibrating  systems. 

An  excellent  treatment  of  the  science  of  friction  and  wear  is  found  in  a 
book  Friction  and  Wear:  Calculation  Methods  (Ref.  1)  by  the  Russian  authors 
Kragelsky,  Dobychin  and  Kombalov.  The  chapter  on  "Dry  and  Boundary  Friction" 
should  be  of  special  interest  to  researchers  in  this  field  and  has  an 
excellent  treatment  of  the  subject  beginning  as  far  back  as  1508  when  Leonardo 
da  Vinci  established  the  simple  law  of  friction  which  assumed  the  friction 
force  to  be  proportional  to  the  normal  load.  Researchers  should  also  find  the 
three  volumes  of  the  book  on  Friction,  Wear,  Lubrication  (Ref.  2)  edited  by 
Kragelsky  and  Alisin.  Another  book,  also  by  Russian  authors  Panovko  and 
Gubanova  entitled  Stability  and  Oscillation  of  Elastic  Systems  (Ref.  3)  treats 
the  subject  of  self-induced  oscillations  with  dry  friction.  In  addition,  the 
proceedings  volume  of  a  NASA-sponsored  symposium  on  Interdisciplinary  Approach 
to  Friction  and  Wear  (Ref.  4)  edited  by  Ku  contains  material  which  cove-s  the 
several  aspects  of  surface  topography,  friction,  and  adhesion,  wear,  etc. 

These  books  along  with  the  papers  that  will  be  referred  to  in  the  sections 
below  have  served  as  the  basis  for  the  viewpoints  presented  in  this  report. 
Because  of  the  interdisciplinary  nature  of  the  subject  matter,  no  claim  will 
be  made  that  all  the  relevant  literature  has  been  surveyed  here.  The  reader 
may,  therefore,  draw  to  the  authors'  attention  relevant  books  or  papers  not 
included  in  this  preliminary  survey  so  that  a  more  complete  list  of  references 
may  emerge  at  the  end  of  this  research  effort. 


2.1  Interface  Characteristics 

Surfaces  of  real  components  used  in  engineering  practice  are  never 
ideally  smooth  so  that  when  they  come  into  contact  with  each  other,  the 
contact  cannot  be  continuous  and  only  certain  parts  of  the  surface  will  carry 
the  applied  load  (see  Fig.  1).  Thus,  the  true  contact  area  is  the  sum  of 
those  parts  of  the  surfaces  where  the  interaction  between  the  bodies  is 
strong.  The  true  contact  area  is  the  result  of  penetration  or  crushing  of 
individual  asperities,  and,  therefore,  the  contact  area  increases  with 
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increasing  deformation.  The  real  contact  area  is  much  smaller  than  the 
apparent  area  of  contact  and  can  be  as  low  as  ".01  to  . 1 %  of  the  apparent 
contact  area"  (see  Ref.  2). 


PE  Ac  AREA 


SECTION  A-A 


Figure  1.  Contact  Areas 


The  research  pertaining  to  the  mechanism  governing  the  formation  of  true 
contact  area  goes  back  to  1939  when  Bowden  and  Tabor  (Ref.  5)  determined,  on 
the  basis  of  measurement  of  electrical  conductivity,  that  the  true  contact 
area  changes  iu  direct  proportion  to  the  normal  load  N .  In  elastic  contact, 
-his  was  shown  to  be  A  »  N^/3.  Further  generalization  for  practical  applica¬ 
tion  in  all  cases  of  contact  between  metallic  surfaces  led  to  the  relationship 
A  =  N/Oy  where  n '  =  yield  stress  of  the  soft  material.  This  law  is  considered 

to  be  applicable  only  when  the  tips  of  spherical  asperities  lie  at  the  same 
leve  l . 

V.  A.  Zhuravlev,  (Ref.  6)  whose  research  work  done  during  World  War  IT 
remained  unknown  for  a  long  time,  calculated  the  true  surface  area  in  elastic 
contact  of  two  rough  bodies.  He  modeled  spherical  segments  of  equal  radius 
for  the  microasperities  whose  tips  were  assumed  to  be  linearly  distributed  and 
obtained 


A  =  K 


(^) 


10/11 


,10/11 
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where  K  is  a  coefficient,  E  is  modulus  of  elasticity,  and  v  is  Poisson  s 
ratio.  Greenwood  and  Williamson  (Ref.  7)  derived  the  following  relationship 
assuming  elastic  contacts. 


where 


A  =  iryA  roFj(h) 


density  of  asperities 
nominal  contact  area 
radius  of  asperities 

standard  deviation  of  height  distribution  of  asperities 
d/  a 


Fn(h)  = 


( s-h)rV>  ( s)ds 


where  9  represents  the  normal  density  distribution  of  heights  of 
asperities,  (i.e.  <f>  is  the  height  distribution  scaled  to  make  its  standard 
deviation  unity),  and  d  represents  the  distance  between  the  reference  plan'..1 
in  each  contacting  surface  (see  Fig.  2).  It  is  assumed  that  all  asperity 
summits  have  the  same  radius  r  and  their  heights  vary  randomly.  Figure  3 
shows  the  manner  in  which  the  contact  area  may  vary  when  the  interfaces  are 
under  "very  high  pressure".  It  may  be  noted  that  the  true  area  of  contact 
under  "moderate  loads"  can  be  less  than  507=  of  the  possible  total  area. 
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Figure  2.  Surface  Characterization 
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Figure  3.  Behavior  of  Contact  Under  Very  High  Pressure 


bden  and  Pires,  in  Ref.  8  also  present  a  simple  argument  to  denonstrat* 
that  the  contact  area  should  vary,  approximately,  linearily  with  the  normal 
load.  In  general,  the  contact  pressure  will  be  a  constant  independent  of  ti 
deformation.  The  total  contact  area,  Ac ,  can  be  represented  as 


A 
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A 


n 


(*4) 


where  AK  is  the  contact  area  at  point  K  and  n  is  the  total  number  of  contact 
points.  Then  the  normal  force  at  point  K  is  given  by 


(5) 


where  pQ  is  the  yield  pressure  for  the  material.  Substituting  for  the  areas 
in  Eq .  ( A)  , 


A 
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(6) 


where  N  is  the  total  normal  force. 


Figure  4  shows  the  nature  of  deformation  of  the  asperities  as  the  normal 
load  increases.  The  symbols  in  Fig.  4  refer  to  the  height  of  the  asperities 
in  microinches  below  which  a  certain  percentage  of  the  surface  lies.  As  the 
load  increases,  the  highest  points  at  the  interface  flatten  as  expected.  But 
the  upward  slopes  in  the  data  show  that  heights  of  asperities  actually 
increase  indicating  an  elastic  deformation  of  the  entire  surface.  Clearly, 
any  law  of  friction  at  interfaces  must  account  for  the  deformation  taking 
place  in  the  neighborhood  of  a  point  in  question. 


Figure  4.  Behavior  of  Aluminum  Surface  Under  High  Pressure  Contact 


Most  of  the  researchers  who  have  tried  to  characterize  the  asperities 
appear  to  agree  on  the  spherical  segment  as  the  most  convenient  geometry  to 
model  surface  roughness.  The  preference  for  this  simple  geometry  is  obvious 
in  view  of  indeterminate  stress  distribution  calculated  for  other  geometries 
examined,  i.e.  a  bar  model  and  a  cone  model  as  shown  in  Fig.  5. 
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Figure  5.  Comparison  of  Pressure  Distribution  Curves  Obtained  in  the  Contacts  of 
Various  Geometric  Shapes 

If  an  ellipsoidal  geometry  is  considered  appropriate,  a  radius  for  an 
equivalent  sphere  is  calculated  from  re^  =  y/ r iongrtran  ' 

2.2  Characteristics  of  Friction  Forces 

A  great  variety  of  phenomena  occur  in  rubbing  contact  between  surfaces  of 
components.  This  is  because  of  the  dual  nature  of  friction  by  which  surface 
friction  overcomes  molecular  interaction  forces  between  surfaces  and 
mechanical  resistances  associated  with  changes  in  the  profile  of  the  surface 
layer.  Kragelsky  et  al  (Ref.  1)  point  to  data  that  indicates  molecular 
attraction  at  an  interface  could  contribute  substantially  to  the  observed 
frictional  forces.  As  can  be  observed  from  Fig.  6,  the  contribution  to 
friction  from  molecular  forces  can  dominate  particularly  for  lighter  load 
conditions.  It  is  assumed  that  a  sublayer  of  decreased  resistance  is  formed 
at  the  interface  between  rubbing  components.  One  speaks  of  dry  friction 
forces  when  the  sublayer  is  in  a  solid  phase. 
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Figure  6.  Variation  of  Friction  Coefficient  Components  with  Normal  Load 


Mechanical  deformation  may  occur  due  to  either  penetration  of  surfaces 
with  ploughing  or  continuous  formation  and  rupture  of  welds.  As  stated 
earlier,  all  these  processes  are  influenced  by  the  environment.  As  observed 
by  Kragelsky,  et  al  (Ref.  1)  "in  spite  of  the  complex  nature  of  mechanical, 
physical,  and  chemical  processes  involved,  it  is  still  possible  to  pinpoint 
certain  principles  which  are  common  to  all  friction  and  wear  processes  .  They 
are  outlined  below: 

1.  Three  interconnected  processes  occur  simultaneously  in  sliding  contact 

a)  Interaction  of  the  surfaces 

b)  Changes  in  the  surface  layers  and  films 

c)  Rupture  of  the  surface  layer. 

2.  The  contact  between  solids  occurs  at  discrete  points  because  of  surface 
roughness . 

3.  Increase  in  load  leads  to  increase  in  contact  area  due  to  increase  in 
number  of  contact  points.  Therefore,  the  real  pressure  at  any  point  increases 
only  slightly. 
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4.  Tangential  resistances  are  additive.  Dissipation  of  energy  in  the 
formation  and  rupture  of  an  individual  frictional  bond  is  determined  by  the 
resistance  in  overcoming  molecular  interaction  at  points  of  contact,  and  the 
mechanical  ploughing  effects. 

5.  Surface  sliding  requires  formation  and  rupture  of  a  thin  layer  (weakened 
layer  of  base  material  or  films)  whose  shear  strength  must  be  less  than  the 
shear  strength  of  the  substrate. 

The  apparent  simplicity  of  the  relationship  that  friction  force  induced 
at  a  point  on  an  interface  is  proportional  to  the  normal  load  at  that  point 
has  been  the  primary  reason  for  its  use  in  most  studies  of  dry  friction. 
However,  in  addition  to  the  problem  of  obtaining  a  reliable  estimate  of  the 
constant  of  proportionality,  i.e.  the  coefficient  of  friction  p,  the  repre¬ 
sentation  F  =  uN  leads  to  computational  complications  even  when  applied  to  the 
study  of  a  single  degtee-of-f reedom  system  (see  Fig.  7),  as  shown  by  Den 


Figure  7.  Den  Hartog’s  Model 


Hartog  (Ref.  9).  Oden  and  Pires  (Ref.  8)  have  recently  discussed  the  mathe¬ 
matical  and  physical  difficulties  associated  with  such  a  representation. 
Nevertheless,  a  substantial  number  of  attempts  have  been  made  and  continue  to 
be  made  to  use  modifications  of  this  basic  model  as  shown  in  Fig.  8  (Refs.  10 
to  26).  Such  mod i f icat ions  are  adaptations  of  the  basic  model  to  the  analysis 
of  particular  components  under  study.  Adding  more  degrees-of-freedom  to  the 
basic  model  leads  to  additional  complexity  (Ref.  10)  which  can  be  avoided  by 
the  use  of  linearized  versions  of  the  model.  As  shown  by  Jacobsen  (Ref.  11), 
such  linearization  can  provide  a  good  representation  of  system  behavior  if  the 
fr.ction  force  is  small  so  that  it  does  not  distort  the  motion  appreciably. 
Subsequent  analyses,  many  of  which  have  been  devised  to  represent  specific 
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Figure  8.  Dry  Friction  Damping  Models 


vibratory  components,  have  generally  relied  upon  the  basic  model  in  its 
linearized  version.  These  developments  have  mainly  involved  more  degrees-of- 
freedom  and  different  configurations  (Refs.  12,  13,  14,  15).  Improved 
friction  models  have  been  tried,  for  example,  making  the  friction  force 
dependent  on  the  vibration  amplitude  (Ref.  16)  and  by  analyzing  the  inter¬ 
facial  slip  in  a  vibrating  beam  (Ref.  17).  However,  there  have  been  few 
studies  which  incorporate  a  detailed  model  of  interface  behavior  into  a  model 
of  a  vibrating  system. 

In  1778,  a  two  term  formula  was  proposed  by  Coulomb  according  to  which 
the  friction  force  F  was  written  as  F  *  AN+B,  where  A  can  be  identified  as  the 
familiar  coefficient  of  friction.  With  such  a  representation,  the  coefficient 
of  friction  is  not  a  property  depending  only  on  the  materials  in  contact  but 
also  on  the  normal  load. 

Interface  behavior  is  a  difficult  area  of  study  which  is  why  it  is  often 
avoided  by  invoking  a  constant  overall  coefficient  of  friction.  Any 
discussion  of  friction  in  the  context  of  vibratory  rub  must  recognize  the  need 
to  distinguish  between  the  tangential  forces  of  external  static  friction  and 
those  of  external  dynamic  friction.  The  former  exist  in  a  region  of  small 
partially  reversible  displacements  and  peak  at  the  boundary  of  such  a  region. 
Dynamic  friction  forces  are  not  dependent  on  the  magnitude  of  displacement. 

The  fully  developed  static  friction  force  is  characterized  by  a  coefficient 
which  is  different  from  that  corresponding  to  the  dynamic  friction  force.  In 
the  solution  of  vibration  problems,  prior  to  the  initiation  of  any  sliding 
motion,  the  static  coefficient  of  friction  should  govern.  The  coefficient  of 
friction  is  either  chosen  on  the  basis  of  practical  experience  or  used  as  s 
disposable  parameter  to  obtain  the  best  agreement  between  theory  and 
experiment.  However,  a  deeper  understanding  of  the  effects  of  interface 
damping  on  vibrating  systems  can  only  be  obtained  by  studying  the  friction 
processes  in  more  detail. 

Introducing  the  concepts  of  a  local  coefficient  of  friction  and  combining 
it  with  a  study  of  stress  distribution  in  the  contact  leads  to  a  microslip 
solution  in  which  only  a  part  of  the  contact  undergoes  relative  motion. 

Goodman  (Ref.  18)  has  reviewed  the  basis  of  this  approach  which  has  since  been 
applied  to  practical  systems  through  the  use  of  finite  element  techniques 
(Refs.  19,  20,  21). 

Postulating  an  overall  friction  force,  but  allowing  it  to  vary  throughout 
each  vibration  cycle  leads  to  the  study  of  stick-slip  behavior.  Antoniou  et 
al  (Ref.  22)  have  reviewed  this  approach  and  show  the  degree  of  complication 
which  this  can  add  to  a  one  degree  of  freedom  system. 


12 


The  classical  mathematical  model  of  dry  friction,  i.e. 


\  +  Fo  ;  *  >  0 

F  =  (6) 

|  -  F0  ;  k  <  0 

(where  Fq  is  the  force  of  slipping  friction,  x  is  the  velocity  of  slipping) 
leads  to  a  force-displacement  relationship  for  alternating  input  displacement 
illustrated  by  two  fixed  horizontal  lines  and  two  vertical  lines  as  shown  in 
Fig.  9.  Experimental  observations,  however,  suggest  that  the  transition  from 


Figure  9. 


one  equation  of  slip 
by  Eq.  (6),  and  more 


(-F0)  to  the  other  (+F0)  is  smoother  than  indicated 
like  the  characteristics  shown  in  Fig.  10  (see  Ref.  27). 


Figure  10.  Hysteresis  Characteristic  of  Dry  Friction 


This  important  observation  allows  the  introduction  of  hysteresis  character¬ 
istics  of  contact  interaction  and  may  serve  as  a  basis  of  the  analysis  of 
friction  damping.  Accordingly,  if  x  and  F  represent  displacement  and  friction 
force  respectively  and  x0,  FQ  their  maximum  values,  the  hysteresis  repre¬ 
sentation  may  be  represented  as 


-BU-Ci) 

1  -  e  1  ;  Z  >  0 

BU-Co) 

-  1  +  e  Z  ;  Z  <  0 


where 


a  =  F/Fq 
Z  ~  x/xQ 
8  =  bxD 

b  *  is  a  material  constant 


constants  Cj,  C2  can  be  shown  to  be  equal  and  opposite  to  each  other. 


i.e. ,  C, 


da  _ 
d? 


8 ( 1-a)  ;  Z  >  0 


8( 1+a)  ;  Z  <  0 


Given  initial  values  of  Z,  a,  the  constants  Cj^  ,  C2  can  be  generated  for  the 
ascending  and  descending  parts  of  the  loop.  The  procedure  can  be  repeated  for 
several  vibration  cycles. 

An  amplitude  dependent  model  of  friction  can  be  derived  as  shown  in 
Ref.  23  for  any  prescribed  value  of  8.  For  example,  for  Z  =  A,  from  Eq.  (7) 


aA  =  1  "  e 


-8(A-C  ) 


a-A  =  +  e 


8 ( -A+C  ) 


which  leads  to 


C  =  ~  In  (cosh  8A) 

°  6 
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using  which  a  can  be  shown  to  be 


a 


1  - 
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-1  + 


cosh  BA 

egg 

cosh  BA 


e  >  o 


s  <  o 


(id 


Equation  (11)  is  an  amplitude  dependent  friction  model.  The  area  of  the 
hysteresis  loop  due  to  friction  can  be  shown  to  be  equal 


4(  A  -  —  tanh  BA) 

B 


(12) 


from  which  6  can  be  determined  from  tests.  It  may  be  noted  that  4A  in  the 
above  expression  is  the  area  of  the  loop  consistent  with  the  friction  charac¬ 
terization  shown  in  Eq.  (6). 

An  appreciation  of  the  continuum  nature  of  bodies  in  contact  leads  to  the 
concept  of  nonlocal  laws  of  friction  postulated  by  Oden  and  Pires  (Ref.  8). 
With  such  a  representation,  the  friction  force  at  a  point  x  is  the  sum  total 
of  the  influence  of  normal  stresses  over  a  neighborhood  as  shown  schematically 
in  Fig.  11.  Thus,  unlike  the  local  law 


Figure  11.  Nonlocal  Nonlinear  Laws  of  Friction  (Oden  &  Pires) 

x(x)  =  yon(x)  (13) 


a  nonlocal  law  represents 


(  14) 


t(x)  =  u  [  Wp(x-y)  an(y)dAc 


where  Wp  is  a  weighting  function  that  distributes  the  influence  of  on  in  a 
circle  of  radius  p.  Recognition  of  the  fact  that  the  friction  (shear  stress) 
vector  is  always  opposed  to  the  direction  of  relative  motion  leads  to  the 
modification  of  Eq.  (14)  as  follows, 


t(x)  =  - un j  /  Wp(x-y)  on(y)dAc 
A„ 


where  u^  is  the  vector  of  relative  displacement  along  the  surface.  A  more 
complete  representation  is  achieved  by  introducing  the  requirement  that  the 
friction  force  be  dependent  on  the  motion  (displacement  or  velocity)  and  leads 


t(x)  =  -Unf<J>e(uj)  /  Wp(x-y)  on(y)dA( 

A„ 


Figures  12  and  13  show  possible  representation  of  W  and  6. 
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Figure  12.  Nonlocal  Effects 
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Figure  13.  Nonlinear  Effects 


2.2.1  Coefficient  of  Friction 


An  accurate  calculation  of  friction  forces  induced  at  interfaces  in 
relative  vibratory  motion  depends  on  the  accuracy  with  which  the  coefficient 
of  sliding  friction  can  be  estimated.  As  observed  earlier,  there  appears  to 
be  no  experimental  data  obtained  in  a  vibratory  environment  from  which  the 
nature  of  the  coefficient  of  friction  can  be  obtained.  The  following  is  a 
summary  of  such  an  experiment  performed  at  the  United  Technologies  Research 
Center  under  contract  to  NASA  LeRC  (Ref.  14). 

The  approach  that  was  chosen  in  that  investigation  was  to  provide 
sinusoidal  relative  motion  between  two  test  pieces  which  were  held  together 
with  a  constant  normal  force  and  to  measure  the  resulting  frictional  forces 
directly.  The  range  of  normal  loads,  frequencies  and  relative  velocities 
imposed  were  consistent  with  vibration  of  a  jet  engine  fan  in  which  resonant 
vibration  is  partly  limited  by  damping  available  due  to  rubbing  at  shroud 
interfaces.  The  test  pieces  were  a  pair  of  titanium  alloy  (8-1-1)  bars  with 
their  rubbing  interfaces  flame  sprayed  with  tungsten  carbide.  The  size  of  tlu 
rubbing  surface  was  22.35  x  6.35mm.  The  test  assembly,  shown  in  Figs.  14  and 
15  uses  an  electrodynamic  shaker  to  induce  sinusoidal  excitation  of  the  lower 
test  piece.  The  upper  piece  which  rubs  on  the  lower  piece  during  test  was 
held  in  such  a  way  that  longitudinal  forces  on  the  piece  could  be  measured 
using  a  strain  gaged  load  bolt,  without  any  appreciable  movement.  The  normal 
load  on  the  joint  was  applied  bv  setting  the  required  weights  on  the  loading 
platform.  Testing  consisted  of  setting  the  required  frequency  and  input 


Figure  15.  Friction  Test  Assembly 


acceleration  and  recording  the  frictional  force-slip  loop  on  an  oscilloscope 
for  a  range  of  normal  loads. 

The  nominal  condition  for  testing  was  defined  as  a  displacement  of 
0.127mm  DA  at  a  frequency  of  280Hz  with  normal  loads  ranging  from  about  250  to 
700N.  Three  loads  were  used,  based  on  the  weights  available  and  the  dead  load 
of  the  upper  arm;  these  were  271,  492,  and  672N.  The  displacement  was  varied 
from  0.056  to  0.132mm  DA  and  frequencies  of  80,  140,  280  and  420Hz  were  input. 
Two  typical  friction  force-slip  loops  are  shown  in  Fig.  16. 

The  equivalent  friction  coefficient  was  determined  from  the  measured  area 
of  the  recorded  (photograph  of  oscilloscope  picture)  loop  and  a  knowledge  of 
the  input  frequency  and  displacement.  The  resulting  values  for  friction 
coefficient  are  shown  plotted  against  normal  load,  frequency  and  velocity,  in 
Fig.  17.  A  least  squares  fit  of  an  equation  of  the  form:  coef  =  (A/Load)+b 
was  attempted  using  all  the  data.  The  resulting  curve  is  shown  superimposed 
on  the  test  data  points  in  Fig.  17.  An  examination  of  the  data  in  Fig.  17 
indicates  that  there  is  no  obvious  correlation  with  either  input  frequency  or 
maximum  velocity. 


2.3  Dry  Friction  Damping  Technology 

Control  of  resonant  vibration  in  light,  flexible  structures  has  been 
found  to  be  largely  the  result  of  damping  available  at  mating  surfaces.  The 
potentials  of  such  damping  depends,  as  has  been  observed  earlier,  on  a  host  of 
parameters  that  include  the  magnitude  of  normal  force  holding  the  surfaces  in 
close  contact  and  the  characteristics  of  the  surfaces.  These  and  other 
parameters  determine  the  nature  and  extent  of  slip  at  the  interfaces  so  that 
an  optimum  in  energy  dissipation  can  be  expected  with  an  appropriate  combina¬ 
tion  of  normal  load  and  relative  motion  (Ref.  23).  Therefore,  an  under¬ 
standing  of  the  processes  involved  in  relative  vibratory  mot  ion  at  joints 
becomes  a  basic  requirement  that  may  influence  the  design  of  structures  where 
friction  becomes  an  important  source  of  damping. 

The  analytical  approaches  for  estimating  friction  damping  can  be  broadly 
classified  as  follows:  (1)  macroslip  approach  and  (2)  microslip  approach.  In 
the  macroslip  approach,  the  entire  interface  is  assumed  to  be  either  slipping 
or  stuck.  The  friction  mechanism  is  either  replaced  by  an  equivalent  linear 
viscous  model  or  assumed  to  be  governed  by  some  form  of  Coulomb's  law  of  dry 
friction.  The  effect  of  friction  damping  is  obtained  by  determining  the 
forced  response  of  the  component.  The  analysis  involved  is  relatively 
straightforward  and  the  justification  for  wide-spread  use  of  this  approach  is 
its  effectiveness  at  predicting  the  actual  response.  However,  there  is  some 
question  as  to  the  validity  of  the  macroslip  approach  when  the  interface  has  a 
large  area,  and  is  subject  to  a  nonuniform  load  distribution.  In  the 
microslip  approach,  a  relatively  detailed  analysis  of  the  stress  distribution 
at  the  interface  is  carried  out,  typically  via  a  finite-element  procedure. 

The  extent  of  local  slip,  not  necessarily  throughout  the  interface,  betwn  n. 
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Figure  16.  Typical  Measured  Friction  Force — Slip  Loops 
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Figure  17.  Variation  of  Equivalent  Friction  Coefficient  During  Sinusoidal  Motion 
With  Normal  Load,  Frequency  and  Maximum  Relative  Velocity 
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pairs  of  contacting  points  is  determined  by  applying  Coulomb's  law  of  friction 
to  the  normal  and  tangential  stresses.  A  detailed  knowledge  of  interface  slip 
dynamics  can  be  obtained  from  this  approach.  Although  much  useful  information 
can  be  gained  in  a  linearized  analysis,  a  fully  nonlinear  time  history 
solution  may  be  necessary  and  feasible  with  modern  computational  methods  and 
equipment . 

Calculation  of  energy  dissipation  over  a  period  of  vibration  requires  an 
accurate  assessment  of  both  the  forces  and  motion  of  an  interface.  Much 
effort  will  be  needed  before  satisfactory  analytical  models  become  available 
to  calculate  these  parameters.  In  the  meantime,  the  process  of  developing  and 
calibrating  analytical  models  requires  advances  in  instrumentation  that  can 
measure  extremely  small  vibratory  motion  at  interfaces.  A  first  step  towards 
such  development  has  been  taken  at  the  United  Technologies  Research  Center  and 
is  described  below. 

The  measurement  system  comprises  a  gallium  arsenide  photoemit t ing  diode 
(Texas  Instruments  Type  TIL24)  and  silicon  position  sensing  detector  (UDT  Type 
PIN-SC/4D)  pair  mounted  side  by  side  in  small  aluminum  block  housings  (see 
Fig.  18).  The  position  reference  is  provided  by  the  infrared  LED  emitting  a 
light  beam,  the  centroid  of  which  is  sensed  by  the  UDT  two-axes  position 
sensing  photo-detector.  Each  axis  of  the  detector  has  two  outputs  which  are 
amplified  to  provide  adequate  signal  level.  The  voltage  difference  between 
the  outputs  of  each  axis  provide  the  X  and  Y  position  information. 
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Figure  18.  Shroud  Motion  Measurement  Device 


A  feasibility  study  of  this  system  was  conducted  to  determine  its 
potential  for  measuring  extremely  small  motions  at  a  shroud  interface.  A 
single  part-span-shrouded  fan  blade  was  used  and  bench  tests  were  conducted 
over  a  range  of  normal  loads  and  frequencies.  The  installation  is  shown  ir, 
Fig.  19. 

The  sensor  assembly  was  glued  using  epoxy  adhesive  onto  the  shroud  of  a 
part  span  shrouded  titanium  fan  blade  which  had  its  root  welded  into  a  massive 
titanium  block.  The  test  item  was  enclosed  in  an  aluminum  fixture.  Applica¬ 
tion  of  the  required  normal  loads  on  the  shrouds  was  accomplished  via  loading 
platens  mounted  on  a  steel  loading  disk.  The  loads  were  set  by  applying  a 
moment  to  the  disk  and  then  locking  the  disk  in  place.  The  diode  assembly  was 
mounted  on  the  load  platen  opposite  the  sensors  with  a  gap  of  12  mils  between 
the  housings.  The  detector  was  oriented  such  that  the  Y-axis  was  radial  and 
the  X-axis  directed  along  the  shroud  interface.  The  whole  assembly  was 
mounted  on  a  slip  plate  attached  to  an  electrodynamic  shaker.  Sinusoidal 
motion  was  applied  to  the  blade  horizontally  in  the  direction  of  flap  of  the 
blade  tip.  The  output  from  the  motion  detectors  was  passed  through  a  low  pass 
3kHz  filter  and  then  displayed  on  a  Tektronix  DM  43  dual  channel  oscilloscope. 
Permanent  records  were  made  by  photographing  the  resulting  wave  forms  and 
Lissajous  figures  for  each  detector  set. 

The  results  of  the  blade  tests  show  some  interesting  features  of  shroud 
motion  for  this  particular  test  arrangement.  The  tests  indicate  that: 

a)  Stick-slip  or  Coulomb  damping  is  not  evident  from  the  traces.  (A 
sample  trace  is  included  as  Fig.  20.) 

b)  Radial  displacement  generally  exceeds  that  along  the  interface. 

c)  The  strain  response  of  the  blade  increases  with  increase  in  shroud 
load  indicating  lower  damping  and  yet  the  relative  motion  between  the  shrouds 
also  increases. 

The  series  of  tests  conducted  on  the  single  fan  blade  demonstrated  that 
the  optical  system  is  very  suitable  for  measuring  small  vibratory  motions.  It 
has  been  established  that  motions  as  small  as  50  x  10-fc  inch  can  be  measured 
with  this  system  in  the  frequency  range  of  d.c.  to  2kHz.  The  results  of 
calibration  indicate  the  system  to  be  adequate  for  the  range  of  frequencies 
and  displacements  that  are  of  practical  interest.  This  instrumentation  system 
has  been  successfully  used  in  measuring  the  extremely  small  vibratory  motions 
at  ten  successive  interfaces  of  a  part-span  shrouded  fan.  It  is  anticipated 
that  in  any  experimental  research  in  this  area,  this  type  of  optical  system 
will  be  found  to  be  invaluable. 
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Figure  20.  Recorded  Shroud  Motions  for  First  Mode  Blade  Vibration 
and  Applied  Shroud  Normal  Load  of  34  lbs 
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3.  RESEARCH  ISSUES  AND  PLANS 


An  accurate  assessment  of  damping  due  to  friction  will  obviously  depend 
upon  the  accuracy  with  which  the  interfaces  and  the  environment  are  charac¬ 
terized.  On  the  basis  of  observations  made  earlier,  five  important  inter¬ 
disciplinary  aspects  may  now  be  outlined. 

1.  Surface  Science:  involving  a  definition  of  the  interacting  surfaces  in 
terms  of  (a)  geometry,  (b)  bulk  material  properties  such  as  moduli  of 
elasticity,  hardness,  etc.,  (c)  surface  topography  and  contact  area,  and  (d) 
surface  chemistry. 

2.  Environment :  involving  (a)  definition  of  mechanical  and  thermal  loading, 
(b)  contaminants,  gases  and  liquids  present,  (c)  humidity,  and  (d)  electrical 
and  magnetic  fields. 

3.  Laws  of  friction  that  relate  1  and  2  above. 

4.  Calculation  of  response  of  a  machine  element  using  1,  2  and  3  above  along 
with  a  knowledge  of  structural  properties  such  as  mass,  stiffness,  etc. 

5.  Calculation  of  energy  dissipation  using  2  and  4  above. 

These  interdependent  and  influential  factors  may  be  represented  in  the 
form  of  a  friction  damping  pyramid  as  shown  in  Fig.  21.  Determination  of 


Figure  21.  Factors  Influencing  Structural  Response  When  Frictional 
Forces  are  Significant 


those  parameters  that  are  most  influential  in  a  given  application  is  an 
important  task,  although  certain  aspects  of  research  in  the  technology  of 
friction  damping  will  be  of  a  basic  nature.  For  example,  the  general  nature 
of  the  coefficient  of  friction  may  be  established  although  its  magnitude  will 
vary  from  one  pair  of  materials  to  another.  Analytical  modeling  for  both 
single  degree  of  freedom  and  multiple  degrees  of  freedom  must  be  developed 
from  first  principles  and  calibrated  on  the  basis  of  carefully  planned  experi¬ 
mental  investigations.  The  recommended  approach  is  shown  in  Fig.  22. 

The  single  degree  of  freedom  model  that  will  be  examined  both  analy¬ 
tically  and  experimentally  is  shown  schematically  in  Fig.  23.  Similarly,  the 
cantilever  beam  of  uniform  cross  section  (see  Fig.  24)  with  a  friction  joint 
at  the  root  serves  as  a  multiple  degree-of-f reedom  system.  The  analytical 
approach  that  will  be  used  in  this  study  is  described  in  Appendix  A  along  with 
typical  numerical  results  obtained  for  the  SDOF  model.  The  aim  of  the  experi¬ 
mental  effort  in  regard  to  both  those  models  will  be  to  develop  a  vibratory 
data  base  that  can  be  used  to  (a)  calibrate  the  analytical  modeling,  (b) 
assess  the  importance  of  various  parameters  and  (c)  quantify  the  extent  of 
damping  available. 

As  noted  earlier,  an  important  aspect  of  friction  analysis  is  the  deter¬ 
mination  of  the  law  of  friction  that  is  valid  and  appropriate  for  application 
to  a  given  vibrating  component.  The  representation  in  the  form  of  a  nonlocal 
law  affords  the  flexibility  needed  in  this  complex  problem  area.  The  repre¬ 
sentation  suggested  by  Oden  (Ref.  8)  was  primarily  intended  for  the  analyses 
of  contact  problems  and  therefore  does  not  include  provisions  for  gross 
sliding  and  vibratory  motion.  Further,  in  such  a  representation,  the  points 
on  an  interface  of  a  vibratory  component  approach  a  zero  relative  velocity, 
and  thus  no  provision  is  made  for  the  velocity  dependence  of  frictional 
forces.  However,  this  concept  provides  a  rational  basis  for  deriving  more 
general  empirical  descriptions  which  are  appropriate  to  a  vibratory  environ- 


A  possible  modification  of  Oden's  nonlocal  friction  law  is  developed 
below  and  will  be  examined  further  during  the  course  of  this  study.  The 
modification  shown  below  incorporates  the  nonlinear  function  <J>  (see  Eq .  (16)) 
into  the  integral,  i.e., 


x ( x ,  t )  =  -v  -r-  J  J 

M  _»  -CO 


/r(x(  t )  ,: 

\  e(xj ) 
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t  <  t. 
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where  x(t^)  =  Xj 
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ANALYSIS  EXPERIMENT 


Figure  22.  Dry  Friction  Program 


Figure  23.  Schematic  of  SDOF  Model  Test  Rig 
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where  aQ,  Ta,  f.Q,  T^,  P  are  material  constants  and  Tsp  is  the  surface 
temperature  at  Xj. 


In  the  characterization  of  surfaces  between  contacting  bodies,  it 
generally  assumed  that  the  asperities  are  random,  but  generally  tooling 


operations  place  grooves  on  the  surfaces  in  a  preferential  direction.  The 
directional  nature  of  these  machining  grooves  will  result  in  a  friction  fore 
that  need  not  align  itself  against  the  direction  of  relative  motion. 
Zmitrowicz,  in  Ref.  28,  accounts  for  this  by  representing  the  coefficient  of 
friction  in  the  form  of  a  2  x  2  matrix.  For  example,  a  frictional  stress  t 
acting  opposite  to  the  relative  surface  velocities  would  be  given  by 
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where 


x,y  are  in  the  surface  coordinates  with  unit  vectors  i,j,  respectively 

z  is  normal  to  the  surface 

vx,Vy  are  the  relative  surface  velocities 

tv,t  are  the  surface  friction  shear  stresses 

A  .  y 

a  is  the  normal  surface  stress,  and 


v  2  +  v  2 
1  x  y 


when  (a)  the  frictional  shear  stress  is  along  the  grooves  for  a  relative 
velocity  along  the  grooves  and  (b)  the  frictional  shear  stress  is  perpen¬ 
dicular  to  the  grooves  for  a  relative  velocity  perpendicular  to  the  grooves, 
then  the  matrix  p  can  be  shown  to  be  symmetric. 
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APPENDIX  A 


SINGLE  DEGREE  OF  FREEDOM  MODEL 


The  objective  of  this  study  was  to  develop  a  mathematical  model  of  a 
simple  system  incorporating  nonlinear  frictional  elements  for  the  purpose  of 
(a)  assisting  in  the  design  of  a  basic  experimental  rig,  (b)  examining  the 
effects  of  different  nonlinear  frictional  models,  and  (c)  proving  a  solution 
method  which  could  be  applied  to  mult i-degree-of- freedom  systems.  The  basic 
model  selected  was  a  single  spring-mass  system  with  a  frictional  force  applied 
directly  to  the  block  and  input  displacements  (sinusoidal  and  step  functions) 
applied  to  the  end  of  the  spring  as  shown  below. 

N 

X=  X0f(  t) 

It  is  required  to  compute  at  any  time  t,  the  motion  x,  of  the  mass  M 
along  the  axis  of  the  spring.  Acting  on  the  mass  are  an  input  force  equal  to 
the  product  of  the  spring  stiffness  and  the  instantaneous  spring  extension 
(X-x)  due  to  an  input  displacement  at  the  end  of  the  spring,  and  a  frictional 
force  F,  which  retards  the  motion. 

The  equation  of  motion  for  such  a  system  when  the  mass  is  moving  is  shown 
below: 


M  d2x/dt  2  =  -K(X-x)  -  sgn ( x )  F  (Al) 

Let  X  =  X0f(^t) 

Then 

M  d2x/dt2  =  -Kx  +  KXQ f (u>t )  -  sgn(x)F  (A2) 


nondiraensionalizing  Eq.  (A2)  by  letting 


q  =  M  w2x/KX0,  t  =  wt ,  s  =  u)n/w,  wn  =  /K/M 
D  =  F/KXq  and  q  =  d2q/dt2 


(A3) 
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leads  to 


q  =  f(t)  -  s2q  -  sgn(q)  D  (A4) 


The  method  of  solution  applied  to  Eq .  (A4)  is  one  developed  by  Houbolt  in 
Ref.  29.  The  approach  is  to  use  difference  equivalents  for  derivatives  to 
develop  a  recurrence  relation  that  permits  step-by-step  calculation  of  the 
response  and  of  the  forces  acting  on  the  mass.  The  main  feature  of  this 
recurrence  approach  is  that  the  generality  and  physical  aspects  of  the  basic 
equilibrium  relations  of  the  problem  are  preserved  without  loss  of  ease  in 
application.  It  is  also  applicable  to  nonlinear  systems  and  MDOF  systems. 

By  the  use  of  difference  equations,  the  differential  equation  (A4)  may  be 
transformed  into  an  equation  that  involves  deflection  ordinates  at  several 
successive  values  of  time.  Considering  a  cubic  curve  passing  through  four 
successive  ordinates  time  c  apart,  as  shown  below, 


t'.e  following  difference  equations  for  the  derivatives  at  the  (n-l)1-^  and 
nth  steps  may  be  obtained. 


=  (2<n 

+  3%-l 

'  6cV2 

+  9n-3)/6£ 

(A5) 

n- 1 

(qn  -  2q 

n-1  +  *^n- 

,2  )/e2 
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+  9c^n-2 
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(A8) 

Substituting  Eq.  (A8)  into  (A4)  gives  the  expression  for  the  nth  step 
value  of  q, 

q  =  (5q  .  -  4q  ?  +  q  ,  +  e2 f ( T )  -  e2 sgn( q )D) / (2  +  e2s2)  (A9) 


.  .....  • 
For  the  case  of  sinusoidal  input,  substituting  the  initial  conditions  qQ  =  qQ 

=  qQ  =  0  at  the  (n-l)1"^  step  into  Eqs.  (A5)  ,  (A6)  ,  and  (A9),  gives  the  expres¬ 
sions  for  the  "first  step"  values  of  q,  q,  q  as 

q^  =  (Sim  -  sgn(q)D)/(6/e2  +  s2) 

qL  =  3qj/E  (A10) 

q i  =  6q j /e 

with  initial  back  steps 

<J-i  =  _c?l 

and  (All) 

q_2  =  -8qL 

For  the  case  of  an  initial  spring  end  displacement,  the  initial  conditions  are 
qQ  =  qQ  =  0  and  qQ  =  1  .  The  "first  Step"  values  for  q,  q  and  q  are  given  by 

qx  =  (2  +  f ( T )  -  D)/(6/e2  +  1) 

c|l  =  3qj/e  -  e/2  (A12) 

9l  =  6q^/e  -  2 

with  initial  "back  step"  values 

9-!  =  e2  -  qL 

and  (A13) 

q_2  =  6e^  -  8qA 

The  above  expressions  have  been  programmed  in  BASIC  for  use  on  an  HP  9816 
computer  with  the  provision  for  detecting  when  the  block  comes  to  rest 
instantaneously  by  determining  when  the  velocity  falls  below  a  small  value. 
This  provision  is  required  in  the  case  of  a  Coulomb  type  friction  model  where 
the  values  of  static  and  dynamic  friction  coefficients  are  unequal.  Based  on 
preliminary  studies,  two  sets  of  results  have  been  selected  for  presentation. 
Figure  A1  shows  the  results  of  an  initial  step  displacement  applied  to  the  end 
of  the  spring.  Decrease  in  the  parameter  D  =  (pN/KX0)  is  equivalent  to  an 
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Figure  A1.  Variation  of  Peak  Response  of  SDOF  Coulomb  Friction  Model  with  Initial 
Force  Ratio  for  Step  Input  in  Displacement 


increase  in  displacement  step  for  constant  frictional  force  (pN)  .  It  is 
assumed  here  that  static  and  dynamic  friction  coefficients  are  equal.  Since 
the  system  is  linear,  the  expected  results  can  be  drawn  ranging  from  a  maximum 
response  of  2  for  no  friction  to  zero  for  the  situation  where  the  frictional 
force  exceeds  the  initial  spring  force.  The  numerical  results  match  the 
expected  results  exactly.  Figure  A2  gives  the  time  histories  of  the  motion  of 
the  mass  (q,  q  and  q)  and  frictional  constraint  forces  acting  on  the  mass  when 
a  sinusoidal  displacement  function  (=sin  u>t)  is  applied  to  the  end  of  the 
spring.  The  parameter  D  has  a  value  of  0.6  when  the  block  is  stuck  and  0.48 
when  the  block  is  sliding.  The  forcing  frequency  is  67%  of  the  natural 
frequency.  The  time  histories  indicate  the  expected  stick-slip  motion  and  are 
in  general  agreement  with  those  obtained  for  a  similar  model  using  limit  cycle 
solut ions . 


39 


-  F/KX0  -  0  0 (STATIC  t.O  4fl(DYNAMIC').S  - 


in  cm  m 

Oil  T- 


m  o  in  o  m  o  in 

O  O  •<-  T-  CM  On 


O)  f~-  in  co  —  ^-roinr^c. 


lOl/ODV-b 


-Ol/Hd 


O)  N  IT) 


co  lT  r-  C- 


C-  lO  co*-^  CO  kT)  N  Cj 


..oi/dsia-b 


t.0in3A-b 


84-3-48-4 


Figure  A2.  Single  Degree  of  Freedom  Coulomb  Friction  Model  Motion  (q,q,q)  and 

Constraint  Force  (FR)  Time  Histories  for  Sinusoidal  Spring  End  Displacment 


APPENDIX  B 


APPLICATION  OF  FINITE  ELEMENT  CODES  TO  DRY  FRICTION  DAMPING 


Five  readily  available  finite  element  codes  were  reviewed  for  their 
ability  to  represent  frictional  forces.  These  included:  ABAQUS  (Ref.  30), 
ADINA  (Ref.  31),  ANSYS  (Ref.  32),  MARC  (Ref.  33)  and  NASTRAN  (Ref.  34).  A 
concept  commonly  used  in  finite  elements  is  to  represent  gap  and  friction 
effects  as  one  dimensional  (truss)  nonlinear  finite  elements.  A  gap  element 
would  be  represented  by  a  spring  connecting  two  nodes  which  has  a  very  low 
stiffness  when  the  nodes  are  separated  and  a  high  stiffness  when  the  nodes  are 
contacted.  Figure  B1  illustrates  the  concept.  In  Fig.  Bl,  the  gap  is  assumed 
to  be  open  initially  with  a  distance  between  nodes  of  6.  The  force  in  the 
spring  is  shown  as  a  function  of  the  relative  displacement  between  the  two 
nodes.  If  the  initial  gap  distance  6  is  negative,  then  the  gap  is  initially 
closed.  In  a  similar  manner,  the  frictional  force  can  also  be  represented. 

If  Coulomb  friction  is  assumed,  the  frictional  force  will  be  opposite  to  the 
relative  displacements,  but  will  be  able  to  take  on  a  maximum  value  given  by 
the  coefficient  of  friction  times  the  normal  load  (the  force  in  the  gap 
element  spring).  Figure  B2  summarizes  the  concept.  Some  computer  codes 
consider  the  frictional  force  to  be  represented  by  a  truss  element  which  is 
e last ic-per feet ly  plastic  with  a  yield  stress  proportional  to  the  gap  force. 

The  ADINA  manual  written  in  1977  did  not  contain  any  mention  of  provi¬ 
sions  for  describing  friction,  although  this  can  be  easily  accomplished  if  not 
already  completed.  ABAQUS  treats  gap  opening  or  closing  and  friction  stick  or 
slip  as  a  boundary  condition.  ANSYS  has  a  similar  treatment.  The  gap  and 
friction  characteristics  are  input  as  elements  but  internally  converted  into 
appropriate  boundary  conditions.  NASTRAN  accounts  for  gap  and  friction 
effects  with  nonlinear  elements.  The  MARC  computer  code  previously  considered 
the  gap  and  friction  characteristics  to  be  represented  by  a  special  nonlinear 
element,  but  this  was  later  abandoned  because  of  inherent  numerical  insta¬ 
bility. 

Currently,  the  MARC  code  uses  a  LaGrange  multiplier  concept  for  the  des¬ 
cription  of  gap  and  friction  effects.  All  of  the  above  representations  suffer 
from  the  fact  that  the  opening  and  closing  of  gaps,  and  frictional  stick  or 
slip  is  highly  nonlinear  and  contributes  a  great  deal  of  numerical  insta¬ 
bility.  Some  of  the  frictional  laws  may  be  adapted  for  incorporating  in  the 
MARC  code  to  determine:  (1)  their  physical  accuracy,  (2)  their  numerical 
stability,  and  (3)  their  ease  of  use  and  implementation. 


NORMAL  FORCE 


Figure  B1  One-Dimensional  Gap  Element 


FRICTION  FORCE 


Figure  B2  Two-Dimensional  Friction  Element 
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